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Abstract 

We demonstrate that if the exponent 7 that measures non-smoothness of the 
velocity field is small then the isotropic zero modes of the scalar's triple cor- 
relation function have the scaling exponents proportional to ^/^. Therefore, 
zero modes are subleading with respect to the forced solution that has normal 

scaling with the exponent 7. 
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Typeset using REVT^ 



INTRODUCTION 



Kraichnan's model of passive scalar advection by white-in-time velocity field has 
became a paradigm within which analytic theory of anomalous scaling in turbulence starts 
to appear 0-lO|. The instrument is the perturbation theory around three limiting cases 
where scalar statistics is Gaussian: i) infinite space dimensionality |]I],|5|,|5|,pi| , ii) extremely 
irregular velocity field which corresponds to a smooth scalar field 0]^,^], and iii) smooth 
velocity field (the Batchelor-Kraichnan limit) P|,P|,|i,TT-n,r7|. The perturbation theory is 
regular in the first two cases and the sets of the exponents obtained agree when the limits 
intersect PJ^. The perturbation theory around the Batchelor-Kraichnan limit is singular 



0], only dipole part of three-point correlation function has been found so far |15[. In this 
paper, motivated by P,p!5[|, we find the isotropic part of the triple correlation function and 
show that the leading term has a normal scaling at small scales. 

It is instructive to discuss first the physics involved to understand the significant differ- 
ence between the first two limits on the one hand and the third limit on the other. Since the 
scalar field at any point is the superposition of fields brought from d directions then it follows 
from a central limit theorem that scalar's statistics approaches Gaussian when space dimen- 
sionality d increases. In the case ii), an irregular velocity field acts like Brownian motion 
so that turbulent diffusion is much like linear diffusion: statistics is Gaussian provided the 
input is Gaussian. What is general in both limits is that Gaussianity is rather uniform over 
the scales, the degree of Gaussianity (say, flatness) is independent of the ratio r/L where 
r is a typical distance in the correlation function and L is an input scale. Quite contrary, 
ln(L/r) is the parameter of Gaussianity in the Batchelor-Kraichnan limit |lT^ , p!6|l so that 
statistics is getting Gaussian at small scales whatever the input statistics. The mechanism 
of Gaussianity is temporal rather than spatial in this case: since the stretching is exponential 
in a smooth velocity field then the cascade time grows logarithmically as the scale decreases. 
This is opposite to what one expects from intermittency and anomalous scaling (anticipated 
beyond the Gaussian limits): the degree of non- Gaussianity has to grow downscales. Already 
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that simple reasoning shows that the way from the Batchelor-Kraichnan hmit towards an 
anomalous scaling at non-smooth velocity field is not to be easy. The formal reason for this 
perturbation theory to be singular is that, at the limit, the many-point correlation functions 
have singularity (smeared by molecular diffusion only) at the coUinear geometry - smooth 
velocity provides for homothetic transformation that does not break coUinearity |]T^. Even 
weak non-smoothness of the velocity smears the singularity i.e. strongly influences the so- 
lution in the narrow region near coUinearity; such a situation calls for a boundary layer 
approach introduced into this problem by Shraiman and Siggia H,0 . 



1. GENERAL RELATIONS 

The triple correlation function of the scalar F(ri,r2,r3) advected by white-in-time ve- 
locity fleld satisfles the closed balance equation [0 

(£ + £rf)F3 = -X3 • (1.1) 

Here, X3(ri,r2,r3) is the triple correlation function of the (non- Gaussian) pumping force. 
Actually it depends on differences = — r^. If \rij\ ^ L (where L is the pumping length) 
then Xs — P3 where P3 is the third-order flux. At growing \rij\ the function Xs tends to zero 
on distances larger than L. The quantity Cd in ( |1.1| ) is the operator of molecular diffusion 
= '^(Vf -|- V2 + V|) {k is the diffusion coefficient), and 

£ = -(1/2) ^/C«^(r,,)Vrvf (1.2) 

is the operator of turbulent diffusion. Here 



/C"^(r) = Dr- 



(^2^°/3_^V) + ^rV^ 
^ ^2-7 



:i.3) 



is the eddy diffusivity related to the velocity pair correlation function. Parameter 7 in ( [1.3|) 
is a measure of velocity non-smoothness, < 7 < 2. 

At 7 = 0, the operator £ is singular for collinear geometry — see ( P?T| ) below. That leads 
to an angular singularity in the correlation functions which is smoothed only by diffusion 



which is therefore relevant at all scales . Contrary, at 7 > the operator C is not singular 
at the collinear geometry and therefore the angular singularity is absent as was pointed out 
in [jl5l. Therefore we can omit the diffusive term Cd in (|1 . 1|) in comparison with C. This is 



possible as long as k ^ •yDr'^. 

In the following we believe xs to be an isotropic function of r^j what dictates the symme- 
try of the solution of (|1.1|) . In this case F3 can be treated as a function of three distances ri2, 



ri3 and only. Then the operator C (|1.2|) can also be rewritten in terms of the separations 



C-^^j:rl-%rl;''-''a, + ... (1.4) 

' l>j 

where dots stand for the terms with cross derivatives dijdki- Since xs = P3 ^'^ ^ij L we 
can easily find a solution of (|L1|) in the region <^ L (cf. Using ( p..4|) we get 

{2-^)PsL^ \^ fn2y' fn^y /rasX^l 



fore 



3D{d-l)d-f 



C 



:i.5) 



L J \L J \L 

which we will call the forced solution. Here C is a constant of the order unity. That solution 
satisfies the equation r^j ^ L but not necessarily matching conditions at r^j ^ L. The 
solution of (|1.1|) that satisfies the condition can be written at rij <^ L as follows: 

F = Ffo,e + , (1.6) 

where Zq is a zero mode of the operator of the turbulent diffusion: CZq = 0. We examine 
here different solutions of the equation CZ = 0. Whether the given mode Z contributes the 
correlation function has to be determined from the matching at Vij ~ L which is beyond the 
scope of our paper. Note that we consider isotropic zero modes while dipole zero modes for 
the anisotropic problem with an imposed mean gradient were treated in . 



It is convenient to introduce instead of rjj the following set of variables 

Xi = — cos6', X2 = — sin6', s = ri2ri3 sin 6^ , (1.7) 

where 9 is the angle between and ria and —00 < a;i < 00, < S2 < cxo, < s < 00. 
Note that the variable s (which is the doubled area of the triangle) is the only dimensional 
parameter among s, Xi, X2- The expressions inverse to (|1.71 ) are 



ri2 = J— , ri3 = J^{xl + xl) 



r23 = J-[{x,-ir + xl], cose= , ■ (1.8) 

V ^2 + Xl 

The operator L and both correlation functions Xa and F3 should be invariant under permu- 
tation of points ri, Y2 and r3, that is under permutations of ri2, and r23. In terms of the 
variable z = xi + 1x2 these transformations can be written as follows 

1 1 . ^ 

1^2: z^l-z\ 2^3: 2^—, 1^3: z^l + , 1.9 

z* z* — 1 

1^2 ^3: z^^—, 1^3^2: z^l--, (1.10) 
1 — z z 

where z* is complex conjugated to z. The doubled area of the triangle s is obviously invariant 

under the transformations. 



2. THE CASE 7 = 



Here, we start with the case 7 = 0. Below we treat the dimensionality d = 2. The results 
can be generalized for an arbitrary dimensionality d. The operator ( |1.2D for = 2 and 7 = 
is rewritten in terms of the variables (|1.7| ) as follows 

(2.1) 



CQ = 2Dxi{di + di). 



Let us stress that derivatives with respect to s are absent in Cq. Then a solution of the 
equation £0-^3 = ~X3 can be written as 



F3{S, Xi,X2) 



1 



SttD j-oo jo 



X2 



In 



{x[ - XiY + {x'2 + 3^2) 
(a;'i - XiY + (x2 - 0:2) 



X3(s,x;,x^), (2.2) 



where we used the explicit expression for the resolvent of Laplacian (cf. |T 



We are interested in the behavior of the correlation functions at r^j ^ L which is not 
sensitive to the particular form of the pumping xz- Thus we can choose any convenient form 
of X3 supplying the convergence of the integral in ( p.2|) . We get the following function 



X3 



l + {r!2 + rfs + rh)/L^' 
5 



(2.3) 



Then the correlation function is given by 



SttD 



Xn 



-oo 



X2 



x'2 + 2sL-2(a;f + x'i + 1 - x[) 



In 



{x[ - XiY + + ^2)^ 
(x'l - + (X2 - X2)" 



(2.4) 



Performing the integration over x[ one obtains 



In 



{X2{1 + t) + wf + {xi - l/2f 

{X2\l - t\ + Wy + (Xi - 1/2)2 



where w = •\/3/4 + x^f^ + L'^X2t/{2s) (cf. H^). The asymptotics of F3 at rjj <^ L is 



(2.5) 



2L) 



In 



2L2 



rf2 + r2 + rig + 2v^, 



The expression ( p.6|) can be rewritten as 

L' ft. {xl + xl)[{x. 



31; ri2ri3r23 



+ const . 



(2.6) 



^In 

QD [(xi - 1/2)2 + + v^/2)2]3 



+ const 



(2.7) 



Here, the first term is the forced solution Fforc which can be obtained from ( |1.5| ) at 7 — 0, 
d = 2, C = 3. The second term in ( [^.7D is the zero mode which has logarithmic 
singularities in the points z = 0, z = 1 and z = 00 that is where one of rjj tends to zero. 
Nevertheless the function F3 has no singularity where one of tends to zero as it is obvious 
from the expression (|2.6|) . Further, as follows from ( p.7|) the zero mode Zq has the linear 
term in the expansion over X2- Remembering (|1.7| ) we conclude that at small angles 6 there 
is the term proportional to \6\ in F3. That is just the angular singularity noted above. 



Repeat again that the singularity is smoothed only by diffusion |17 



We see that the zero mode at 7 = does not depend on the dimensional parameter 
s. Besides, any function of s is a zero mode of the operator Cq since it does not contain the 
derivative over s what is a remarkable property of the case d = 2 and 7 = 0. Nevertheless 
all the zero modes do not contribute to Zq. 



3. ASYMPTOTICS AT SMALL ANGLES 



Because of the scaling properties of C, it is possible to seek the zero mode in the form 
Z (X s^. It is convenient to introduce the following function 
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Z=[^^y{l + {xl + xl)^ + [(1 - x,f + xlf} X{xu X2) , (3.1) 



The function X{xi,X2) should be invariant under all transformations ( |1.9| , |1.1CI| ) and have 



no angular singularities since the proportionality coefficient between Z and X is equal to 
^12 + ^li + as follows from ( pTSQ . 

The equation CZ = can be rewritten as CxX = where Cx is a differential operator 
of the second order over di = d/dxi and 82 = 8/8x2- Coefficients at the derivatives are the 
functions of Xi, X2, A which can be found from (|1.1|J1.8| ). The functions are quite complicated. 
Fortunately only particular parts of the operator Cx will be needed for us. 



At 7 = the operator Cx is determined by (|2.1| ). The operator tends to zero at X2 — ^ 0. 
Therefore at small X2 besides we should take into account also the residue. The term 
leading at small X2 can be written as 

£x oc £2 = [2x1 + co{xi)]8l - AAX282 + 2A(A + 1) , (3.2) 

37 

Co(x) = — ^(1 — x)x[x In + (1 — x) In |1 — x|] . (3.3) 

Note that Cq > 0. The expression (|37^ ) is correct if X2 ^ |xi|, |xi — 1|; |xi — 1| 3> 
exp(— 1/7); |a;i|, \xi — 1| <^ exp(l/7). The asymptotic behavior of a solution of C2X = at 



X2 <^ ^/co IS 

X = Ai{xi) + A2{xi)x2 , (3.4) 

where Ai{x) and A2{x) are arbitrary functions. The analyticity of X at small angles excludes 
the second term in ( |3.4| ) since it supplies the contribution to X which behaves oc The 
equation C2X = can be solved explicitly, a solution having the asymptotics ( ^I^) with 
= is 

^ ^ V 2 ' 2'2' co(xi)y ^ ^ 

Here, F{a, jS;'-/] z) is the hypergeometric function. 

The expression ( |3.5D gives the behavior of the zero mode in the vicinity of the boundary 
layer X2 ~ \/co- If we are interested in the behavior of the zero mode outside the boundary 



layer then it is more convenient to return to Z since C can be approximated as (|2.1| ) and 
consequently Z is a harmonic function there. The asymptotics of ( p.5|) valid at 3> 
gives 

/ 2 

Zoc(A + l)cos(7rA/2)-X2sin(7rA/2) . (3.6) 

\co(a;i)/ 



The behavior (|3.6| ) occurs outside the boundary layer but at small X2- Since we are interested 
in small A we can suggest A ^ 1. Thus we come to the following problem: find the harmonic 
function Z{x\, X2) in the upper half-plane 0:2 > at the boundary condition 



^^-^2^ + ^ = 0, (3.7) 

which is imposed on the function Z at X2 = since at small 7 the width of the boundary 
layer is negligible. 

Let us show that there is no zero mode with A ^ y^. A harmonic function Z{xi,X2) 
inside the region can be presented as an integral of its normal derivative dZ/ dn along the 
contour which is the boundary of the region: 

Z{z) = -i\dt\^^\n\z-t\, (3.8) 



TT 



dn 



where z = xi +1x2 and t is the complex variable going along the contour. Let us consider the 
contour consisting of the semi-circles going around the singular points 0, 1, cxd and the parts 
near the real axis (outside the boundary layer but at small X2) that link the semicircles. If 
A ^ then the condition (|3.7|) tells us that in this case only contributions to (|3.8| ) from the 



semicirles are relevant since the contributions from the parts of the real axis are negligible. 
Separate consideration of the vicinities of the singular points 0, 1, 00 (see below) shows that 
the logarithmic derivative of Z has to be bounded there. Thus the only possible contribution 



to the zero mode associated, say, with the singular point z = is oc Relnz = In + 0:2. 
The complete zero mode should be symmetric under transformations (|1.9|). Performing all 



the transformations to In yxf + and summing the results we obtain zero. That means 
that the function possessing the required symmetry does not exist. 
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Another (equivalent) way of showing that there is no zero mode with A ^ is to 
continue Z to negative X2 by Z[x\^X2) = Z{xi, —X2). Becase in our case we can beheve 
d2Z{xi,X2) = the function should be harmonic in semicircles surrounding the singular 
points 0, 1, oo. Then one can use the properties of the analytical functions in the circles to 
exclude the existence of zero modes with bounded logarithmic derivatives near the singular 



points. Note the difference with the dipole case where such mode has been found 15 



4. VICINITIES OF SINGULAR POINTS 



Here, we describe the set of zero modes that do not have additional smallness of A relative 



to so that the whole boundary condition (|3.7| ) is to be accounted. Necessary information 
about the structure of the modes can be extracted from the analysis of the vicinities of the 
singular points z = 0, z = 1 and z = oo where one needs a separate consideration. Using 
the symmetry properties ( |1.9| , |1.10D we can reduce the consideration to the vicinity of one of 
the points, say z = 1. At 0:2 ^ 1 and |xi — 1| ^1 the operator Cx acquires the following 
form 



lx = fi [p'^dl + 3pd, + 3dl] + 2 sin^ [p'^d^ + pdp + 3% 

-4A[sin2 ippdp + cos sin ipd^\ + 2A(A + 1) , 

X\ — \=p cos V9 , X2 = p sin , ~ 2 ^^^^ ■ 



(4.1) 
(4.2) 



In the exponentially narrow vicinity of the singular point where p -C exp(— 1/7) one has 
/X ^ 1 and the equation = is reduced to 



p'dl^'ipdp^'idl 



X = 0. 



(4.3) 



Solutions of the equation ( [4.3| ) can be expanded into the Fourier series over the angular 
harmonics 



Xm oc sin(mv9)p'* 



(4.4) 



Substituting ( [4. 41) into (|4.3| ) we obtain 



A„ = -1 + Vl + 3m2, (4.5) 

where only nonnegative Am are taken since X should remain finite at p ^ 0. Thus besides 
a constant corresponding to m = we obtain the next term corresponding to m = 1 
that behaves oc p. Since inside the exponentially narrow vicinity a zero mode X is a linear 



combination of ( [4.4|) we come to the conclusion that the matching condition on the boundary 
of the vicinity should be imposed on the logarithmic derivative of X (or Z) which remains 
constant there. 

Now, let us consider the region 1 ^ p ^ exp(— 1/7) where 

p=i7ln-<l. (4.6) 
/ p 

Here, we can consider separately small angles (p where CxX = is reduced to 

(3p + 2^^)dl - 4A^9^ + 2A(A + 1)] X = , (4.7) 

what exactly corresponds to C2X = 0. Again an appropriate solution of ( |4.7| ) is as follows 

+ A A 1 2(^^ 



/ 1 + A A 1 2w2\ 
XocF ■ ■-^] , (4.8) 



that gives the asymptotics 



ttA , , 

Zocl--=yp, (4.9) 

A/op 



at 1 ^ ^9 ^ yffl. Substituting now (|4.6| ) we come to the boundary condition 

d^Z = , Z, (4.10) 

/37ln(l/p) 



imposed on the harmonic function Z at if = 0. Note that ( |4.10|) is nothing but the limit of 
( p.7| ) at p ^ 1. That boundary condition is simple enough and permits explicit expression 
for Z near the singularity. 

Let us represent Z in the following form 

Z = Rei^f (^n^ + tif^ + f (\n^ + tTT - icp^ | (4.11) 
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which is obviously harmonic and invariant under — > Ti—ip. Taking into account ln(l/p) ^ 1 
we obtain from ( [4 .111 ) 



Z(p,<^ = 0)=2Re|/(^ln-j| 

|^(p,^ = 0) = 7rRe|r (ini 

Thus we see that /(x) satisfies the following equation 

2A 
a/37x 

Asymptotic behavior of the solution at x ^ 1 is 



(4.12) 
(4.13) 



(4.14) 



exp 



4 /4A 



2\ 1/4 



"3 V 37 



X 



3/4 



(4.15) 



Expanding ( [4.11 ) with ( [4.15 ) over In (1/p) we obtain 



Z oc cos 



4 /4A 



2\ 1/4 



3/4 



In- 



+ 



(4.16) 



3 \ 37 y y p^ 

where 0o is some phase. We can believe |0o| < ^! its actual value has to be determined by 
the matching at — In p ~ 7. 

Symmetry requirement with respect to Xi — 1 — Xi leads to the condition diZ{l/2, X2) = 
which can be used as the quantization rule for the zero modes having the asymptotics 
( [l.lti| ). They can be classified in accordance with the number of zeros n which the function 
Z have where Xi goes from 1/2 to 1 — exp(— 7). Using the expression ([4.16[) we conclude 
that Amin = (^\fl and for n ^ 1 



(4.17) 



with yet unknown numerical factors which are of order unity. Note that nonsymmetric 
zero modes Z (with another values of a,/3) may exist yet they cannot contribute to Zq. For 
all the modes, the dependence A (7) obtained here has an infinite slope at zero which has 
been also observed in numerics |]TB|. We conclude that the set of zero modes thus found at 
small 7 has exponents larger than the exponent 7 of forced solution. Therefore, the isotropic 
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part of the triple correlation function is shown here to have a normal scaling for sufficiently 
small 7. Since at 7 = 2 the lowest zero mode has A = 4, it is likely that the scahng of the 
isotropic part of the triple correlation function is normal for all 7. 
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